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{Nl From the point of view of tunneling, the physical meaning of the Davies critical 

point of a second order phase transition in the black hole thermodynamics is 

<N 

clarified. At the critical point, the nonthermal contribution vanishes so that the 
black hole radiation is entirely thermal. It separates two phases: one with radiation 

00 

enhanced by the nonthermal contribution, the other suppressed by the nonthermal 
J contribution. We show this in both charged and rotating black holes. The phase 

Q_i transition is also analyzed in the cases in which emissions of charges and angular 

momenta are incorporated. 
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I. INTRODUCTION 

Upon the discovery of the thermodynamic interpretation of black holes [1], the first critical 
phenomenon identified was the second order phase transition discovered by Davies more than 
thirty years ago [2]. The phase structure was analyzed by Hut [3], and it has been generalized 
for varying charges of isolated charged black holes[4]. The issue of stability in the Kerr black 
hole case has been addressed[5]. Yet, all these are still in the thermal radiation context based 
on the Hawking process. Here we present a new proposal to clarify the physical meaning of 
this critical phenomenon in the context of the nonthermal radiation of Parikh-Wilczek[6][7]. 
The existence of nonthermal radiation from black holes have been anticipated [8] [9] beyond 
Hawking's thermal radiation[10]. 

Black holes are unusual thermodynamic systems because we cannot see any structural 
changes directly as far as critical phenomena are concerned. The entropy changes smoothly 
so that checking the behavior of the horizon cannot tell us either. So the only way to 
tell what is going on in the black hole phases is if there is any difference in emissions of 
information from different phases. So the tunneling argument of Parikh-Wilczek is a good 
candidate to investigate if there are such emissions telling about the black holes phases. In 
this paper, we find that indeed this is the case, and that, using this information, we can 
explain the physical meaning of the Davies critical phenomenon. 

In short, we find that there is a competition between thermal part of radiation and 
nonthermal part. At the critical point, the nonthermal part vanishes, leaving only thermal 
radiation that peaks. The critical point separates two phases: In one phase, the nonthermal 
contribution enhances the total radiation, while in the other phase the nonthermal contri- 
bution actually suppresses the total radiation. Once we extend to the cases of emissions of 
charges or angular momenta, we can observe over all enhancement or suppression, but the 
peak of emission remains at (or near) the critical point and separation of two phases persists. 
We could check the characteristics of each phase of black holes with respect to the radiation 
and emissions of charges and angular momenta by introducing effective free energies. 

This paper is organized as following. In section 2, we review and identify what the Davies 
critical point is. In section 3, we explain this critical phenomenon in the RN (Reissner- 
Nordtrom) black hole case, using the tunneling argument. Then, in section 4, it is extended 
to incorporate the emissions of charges and angular momenta in the KN (Kerr-Newman) 
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case. Finally, in section 5, some further comments are given. 

II. DAVIES CRITICAL POINT 

The Davies critical point is identified by the singular behavior of specific heat at some 
nontrivial value of charge-to-mass ratio, Q/M, or, angular-momentum-to-mass ratio, J/M 2 , 
away from the extremal limit [2]. A second order phase transition occurs at this critical point 
and the phenomenon is generic for any charged or rotating black holes. The nature of this 
critical phenomenon is not entirely clear except it has been known that the specific heat 
changes the sign abruptly. 

To illustrate the Davies critical point, we shall start with the RN black hole. Fig.l shows 
the relationship between the Bekenstein-Hawking entropy and the Hawking temperature for 
some value of charge Q, and, in this case, the Davies critical point is the turning point 
marked by II. At the critical point, \Q\/M C = y/3/2, T c = 1/(9ttM c ) and S c = (9/4)ttM c 2 for 
a given charge Q. As Q varies, the critical points trace along S c oc 1/T C 2 . As one can see in 
the plot, specific heat, which is related to the slope of the curve, dS/dT, is singular at II, 
indicating the occurrence of some kind of a second order phase transition. The specific heat 
is negative in region I and positive in region III. The critical exponent for the specific heat 
is 1/2 such that, as T# approaches T c , 

\c Q \~{T c -T H )-V\ (1) 

The Helmholtz free energy of RN black hole (fig. 2) shows a cusp formed at the Davies 
critical point. Note that, unlike normal thermodynamic systems, the part with negative 
specific heat of Helmholtz free energy is convex. The Gibbs free energy also displays a 
turning point at the Davies critical point. In the extremal limit, Gibbs free energy vanishes 
for RN black hole (fig. 3), but in the KN case Gibbs free energy no longer vanishes due to 
the rotational degrees of freedom (fig.4). 

Compared with the Schwarzschild black hole case, in which the specific heat is always 
negative and temperature does not have a turning point, in the RN or KN case there must be 
something happening to generate such a turning point that the behavior of the temperature 
changes, at the same time, changing the specific heat from negative to positive. 
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FIG. 1. S vs. T plot of RN black hole 
for fixed Q = 0.1: The turning point II 
is Davies critical point 
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FIG. 2. Helmholtz free energy of RN 
black hole for fixed Q = 1: The end 
point of the cusp is Davies critical point 
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FIG. 3. RN case: At T = 0, G = 0. 



FIG. 4. KN case: At T = 0, G > 0. 



III. DAVIES CRITICAL POINT FROM TUNNELING IN 
REISSNER-NORDSTROM BLACK HOLE 



The appearance of such a charge-to-mass ratio as a critical point away from the extremal 
limit is quite intriguing and it certainly calls for deeper understanding of its physical 
implications. In this paper, We argue that what happens is due to the nonthermal black 
hole radiation. This can be seen by using Parikh-Wilczek's tunneling argument [6] of black 
hole radiation as follows. 
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In tunneling, the emission rate is given by 

T ~ exp(-2Im7) 

= exp ^—2 Im J p r dr 

rM 



-21m J \f3 M (H)dH \ , 



where 

r out ^ r 



(2) 



f3 M (H) = 2 / - (3) 

is related to the inverse Hawking temperature as we shall see later. After integrating, which 
is done by Parikh-Wilczek, we get 

T ~ e- Arff(w)w , (4) 



where 



^-Mu) ee — ( (M a;) 2 + (M — ui)y/(M - u) 2 - Q 2 
In uj \ 

— M 2 — M^M 2 - Q 2 ^j . 



(5) 



T c s = is called the effective temperature, which is interpreted as the indication how 

temperature changes after the emission of radiation. Note that this effective temperature is 
not the temperature after the emission of radiation, but it is a quantity still defined at the 
temperature before the emission. Assuming 1 3> u/M and I 3> (2ooM — oo 2 ) / (M 2 — Q 2 ), we 
can expand /3 c ff(w) as 

^Mco) = — j^- - Q) + h.o., (6) 



where 



T H = = ^ , 2 (7) 

^ (m + v/M 2 - Q 2 J 

is the Hawking temperature and 

3 M 1 M 3 

/(M, Q) ee 1 + - — _ ^ 2)1/2 - - — _ ^ 2)3/2 . (8) 

The key observation here is that the sign of f(M, Q) depends on M and Q. In fact, 
f(M ) Q) = for \Q\/M C = y/3/2, which is nothing but the Davies critical value! Note 
that f(M, Q) ^ for M ^ M c . 
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If we can relate f(M, Q) to the specific heat, we should be able to understand this better. 
Indeed, in terms of specific heat 

1 fdM\ 
CQ =M\Jf) Q 

n (M+^M 2 -Q 2 ) 2 ^M 2 -Q 2 

= 2,71 = 

M 2 - 2M^M 2 - Q 2 



4ttM f(M, Q) ' 



(9) 



we can reorganize eq.(6) as 



M») = Ph+\^ + 0{u?). (10) 



This clearly shows how Davies critical point appears in the tunneling context. As we shall 
show later, this is generic for any black hole with the corresponding specific heat. Now one 
can see that, for cq < 0, /3 c s < (5h so that (effective) temperature rises, and, at the Davies 
critical point, /3 e fr = Ph, but as soon as black hole crosses over the Davies point into the 
region of cq > 0, the temperature turns lower. 

This also means the following: Ahead of the Davies point (from the right side in fig. 5), 
cq < and the (emission rate of) radiation is enhanced by nonthermal radiation on top of 
the thermal radiation. At the Davies point, the specific heat is singular and the radiation 
is purely thermal. But as the black hole crosses the Davies point, cq > and the radiation 
is reduced by the nonthermal part, lowering the temperature. Even though the assumption 
for this tunneling, (M — u) > Q, excludes the extremal case (the solid line in fig.5 starts 
off away from the bottom), but the final formula still seems to make sense in the extremal 
limit, where the emission rate does become zero as we expect. So, we can safely extrapolate 
this to the extremal limit with caution. Note that the thermal radiation peaks at the Davies 
point and the total emission rate peaks at M — M c + uj. Since u is small enough compared 
to M, it is safe to say that the total radiation peaks at the Davies point. 

Compared with the Schwarzschild case, in which there is no such a turning critical point, 
in the charged case, clearly the difference is caused by the charge Q. So the presence of 
charge is the cause behind such a phase transition with nonthermal radiation. Since angular 
momentum behaves similarly to charge for black holes, we can anticipate this structure will 
persist when J ^ 0. 
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FIG. 5. Plots of (normalized) emission rate vs. mass for fixed Q. They are 
exaggerated to visualize the difference clearly. The solid line is for RN black 
hole with nonthermal contribution, while the dashed line is just for thermal 
part of RN black hole. The dotted line is the thermal part of Schwarzschild 
black hole for comparison. As Q increases, the emission rate decreases and the 
peak location moves to the lower-right direction. 



For the purpose of generalization, let's derive eq.(10) more formally. This can be done as 
follows: First, note that in terms of Wick rotation t E = it 



where (3 is the inverse Hawking temperature, which can be shown by directly integrating 
along a geodesic that crosses the horizon. (See the Appendix A for the proof.) The factor 2 
is present because one needs to cross the horizon back and forth to complete the period for 
the Euclidean time to be a temperature. Using this, we can identify 




dt E = -i/3, 



(11) 




(12) 
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Then, expanding in terms of u', from eqs.(2)(4) we can obtain 

/3 cS (co)co= / 0(M - J) dJ 
Jo 



u'=0 



= u/3(M) + \u?"-f +0(u 3 ) 

=^(M) + r 2 'fi M » 2 +0 (^), (is) 



where /3(M) = (3h- Note that in the above derivation we have not specified what kind of 
black hole we use. So, we have a general theorem that for any black hole (including the case 
with J 7^ 0, i.e. the KN black hole) which emits energy u <C M, eq.(10) holds true. 

The Helmholtz free energy cannot be used to check the stability of the system with respect 
to the radiation because it is constant for fixed Hawking temperature for which both Q and 
M need to be fixed. However, we can take an analogy of the idea of effective temperature 
and check the stability with respect to the radiation. For this purpose, let us define an 
effective Helmholtz free energy after the emission such that 

F eS := M cS - T cS S eS . (14) 

This can be used to compare the Helmholtz free energies before and after emission of 
radiation at fixed temperature. Then using the values after emission 

F cS = M + AM--^-(S + AS) = F + AF, (15) 

Peff 

where 

AM = -u < 0, (16a) 
AS = -f3 cS u < 0, (16b) 

af 4^+<^ 2 ) ( i6c ) 

ZCqM 

and F is the usual Helmholtz free energy, hence AF measures the difference at the same 
Hawking temperature. In fact, 

F cS = M - T cS S. (17) 

Note that the change of the Helmholtz free energy is given in terms of the specific heat in the 
leading order of lu/M. In particular, AF = at the critical point. AF < for c Q < 0, which 
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can be interpreted as that the black hole with nonthermal radiation is more stable so that 
black hole will keep emitting nonthermal radiation until it gets evaporated. However, for 
cq > 0, AF > 0, which means the black hole before the nonthermal radiation is more stable, 
hence the nonthermal radiation will be suppressed, or there has to be another mechanism 
for more nonthermal emission. This will be addressed in section 5 again. 

IV. DAVIES CRITICAL POINT FROM TUNNELING IN KERR-NEWMAN 

BLACK HOLE 

Next, we shall consider a more general case with J ^ 0. At the same time we shall allow 
charged particles[ll][12] (also see [13] for an earlier attempt) as well as angular momenta to 
be emitted. The thermodynamic relation of a KN blackhole is given by 

dM = T H dS + V H dQ + n H dJ, (18) 

where 

Th = T = ! - ± f Iz ' (19a) 

S = l -A Kl (19b) 
Vh - (19c) 

A K 

n H = — , (i9d) 

A K 

where a = J/M, A K (M, J, Q) = An(r 2 + + a 2 ) and r ± = M± a/M 2 - Q 2 - a 2 . 

Then, as shown in Appendix B, for emission of energy u, charge q, and angular momentum 
j, we have 

-2ImJ =AS 

= - p H (u - qV H - j^l h) 

2 ( Ph Pi \ 
2 Mc Q 2Q \Qk q q>qQ J 

_ i -2 f Ph_ _ Ph n \ 
2 \Jkj ajJ J 

qPh .(Ph v PI 

-u^ + h.o., (20) 
J a j 
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where the following definitions of thermodynamic quantities for black holes are used: 

aq = ^tt^ : charge expansion coefficient, (21a) 

\a£ (J 1 

1 dJ 

a j = -j-^f '■ angular momentum expansion coefficient, (21b) 
kq = -j- ( T^y-) '■ charge "compressibility", (21c) 

Q \ (jV h J t 

kj = — ( ] : angular momentum "compressibility". (21d) 



Note that 



{2H yH = 2Jl + ' +g ' )>0. (22a) 



Qk q a Q Q V 51/2 M 2 fi 3 / 2 

& 9 l<oifQ>o, 

" S ( 22C ) 



MW2 (>oifg<o, 

#2, /3i Oa > if Q > 0, 

where B = 1 — (Q 2 + a 2 )/M 2 and J > is assumed without loss of generality. This shows 

that the only term which changes sign for nontrivial values of Q and J is the specific heat 

term. In this case the sign change occurs at 

Q 2 1 / a 2 a 4 \ 

'3-6— -— , (23) 



M 2 4 V M 2 .U ; , 
which identifies the Davies critical points [2]. As pointed out by Davies, this critical phe- 
nomenon takes place even for Q = at a 2 /M 2 = J 2 /M 4 = 2a/3 — 3. In between, as a 2 /M 2 
varies from zero to 2\^3 — 3, Q 2 /M 2 varies from 3/4 to zero. 

One can also check the signs of other thermodynamic quantities and find 

a Q < 0, (24a) 

aj < 0, (24b) 

k q > if Q > 0, < if Q < 0, (24c) 

k.j > 0. (24d) 
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FIG. 6. Plots of (normalized) emission rate vs. mass for fixed Q and J. Again, 
they are exaggerated to visualize the difference clearly. The lower solid line is 
for KN black hole with nonthermal contribution for q = = j, and the upper 
solid line is for KN black hole with nonthermal contribution for q ^ ^ j, 
qQ > 0, while the dashed line is just for thermal part of KN black hole. The 
dotted line is the thermal part of Schwarzschild black hole for comparison. As 
Q or J increases, the emission rate decreases and the peak location moves to 
the lower-right direction. 



In the extremal limit, they all approach to zero, including the specific heat, and, as M — > oo, 
they all diverge. So, the existence of the Davies critical point is quite unique for the specific 
heat. 

Even though u,q,j, are introduced independently, they cannot be arbitrary. This is 
because physically charges or angular momenta cannot be emitted without changing energy. 
For example, it is necessary that 



for the emission rate to make sense and that AS < 0. Otherwise, the emission rate diverges 
as temperature goes to zero. This provides us interesting bounds. 
Let's consider j = case, first. Then 




(25) 



q 



> \v H \ = 




(26) 



11 



for qQ > 0. \Vh\ < 1 and approaches zero as M — > oo, hence this does not necessarily 
imply oj > q. Nevertheless, this bound is quite significant, compared with m e /e ~ 1CT 21 or 
rrip/e ~ 1CT 18 . This means that charge has to be emitted with a sufficiently large energy. Not 
surprisingly, the larger Q/M ratio is, the more energetically charges will be emitted. This 
also has a rather interesting alternative interpretation. Suppose Q/M violates this bound 
for some charge emission, then the bad behavior of emission rate can be interpreted as black 
holes inability to sustain charge-to-mass ratio Q/M. There are two ways of lowering Q/M: 
either emit charges or absorb more mass to increase M. What if not much matter is around 
the black hole to absorb? The only resolution will be to emit charges. In other words, this is 
a bound for a black hole how much charge it can sustain. Based on the classical argument, 
we know Q/M < m/q is a limit for a black hole to absorb stationary charged particle. Here 
we have another bound beyond which a black hole cannot emit charges normally, but must 
get rid of them drastically. Note that, unlike the classical bound, the bound in the tunneling 
case is given in terms of emitting energy, not rest mass. If qQ < 0, eq.(25) is always satisfied. 
This means that opposite charges can also be emitted and Q can actually increase, although 
the emission rate is lower than that of qQ > 0. 

A similar bound can be obtained for the emission of the angular momenta. For q = 0, 
since jJ > 0, 

->n H = — — . 27 

J M 3 (l + vAB)2 + J2/ M 4 

To look into the next order contribution, we first need to impose eq.(25). Now the 
emission rate actually depends on the signs of Q and q. The q 2 and j 2 terms always enhance 
the emission rate, but the contribution of u 2 term depends on the sign of specific heat and 
vanishes at the Davies point, where the emission rate peaks. If qQ > 0, the last two terms 
with q in the quadratic expansion enhance the emission rate. However, if qQ < 0, they 
suppress the emission rate. The term with j always enhances the emission rate since we 
assume J > and j > 0. 

In the RN case, we checked the stability of the black hole with respect to the radiation 
of energy, using the effective Helmholtz free energy. In the KN case, however, the effective 
Helmholtz free energy is not sufficient because of charge and angular momentum emissions. 
Note that the Gibbs free energy cannot be used to check the stability of a KN black hole 
with respect to the radiation either because it is constant for fixed T H , Vh and f2#, for which 
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all M, Q, J are fixed. 

So, we need to consider the effective Gibbs free energy analogously defined, like the 
effective Helmholtz free energy, in terms of effective quantities as 

G e fi = M e g — T eS S e s ~ V e f[Q e ff — f2 eff J e ff , (28) 

then the difference between the effective Gibbs free energy and the usual Gibbs free energy 
at the same Hawking temperature reads, as derived in Appendix C, 

AG = G cff - G = tuGu + qG q + jGj, (29) 



where 



7T 



2 



A ( , - 2M 2 \ 

I Q fl/f2 i /02 i \ 



G ^ -pr-m \ 2M + Q+ WT 2 )' (30a) 
ttQ 



G n = 



>H 



B 3/2 



,„ 2tt Q 2 + 2J 2 /M 2 ( , - 2M 2 \ . 
x ( I + S 3 / 2 - — « M3 7 (^2M 2 + Q 2 + ^ ) ) . (30b) 



ttJ / 1 2n2M 2 -Q 2 f 9 - 2M 2 \\ , x 

where £ ee 1 - (Q 2 + J 2 /M 2 )/M 2 as before and 

72 (Q2 t2(Q2 

A = 25 3 / 2 + 2 - 6^- - 3^= + (31) 
M 4 M 2 M 6 v ' 

Note that A = identifies the Davies critical points. One can show that Gj < for J 7^ 
always, however, the signs of G w and G q depend on Q, J and M. 

In summary, there are over all three different cases of AG depending on the variables, 
which are M, Q, J and u, q,j: AG < for all M, AG ^ for M M , where M is another 
solution to AG = 0, or AG > for some interval of M while AG < for the rest M. If 
J = 0, M = M c , but if J 7^ 0, M is different from M c . The existence of opposite signs of 
AG indicates there are two different phases of KN black holes. The details are as follows. 

If J = 0, this describes the RN black hole with charge emission for q 7^ or without 
charge emission for q = 0. In the q = case, AG correctly reduces to AF so that it simply 
reproduces the result of effective Helmholtz free energy case. So we shall consider q 7^ case 
here. For J = 0, G q can be rewritten as 

r _ *Q 2M 2 {l + B" 2 )-Q 2 

q /3 H B 3 / 2 (2M 2 -Q 2 + 2M 2 B 1 / 2 )(1 + B 1 / 2 ) 1 ' 
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so that G q = if A = 0. G w = = G q , hence AG = at M = M c , that is, the Davies 
critical point. For J = and fixed Q, G w ^ for M ^ M c , while G q ^ for for M ^ M c . 
So, if cu/g > 1, AG ^ for M M c . In the case u/q < 1, there is a value g' such that, if 
u/q < u/q' < 1, J = 0, AG > for M c > M > M , while, if u/q' < u/q < 1, J = 0, AG > 
for Mo > M > M c , where Mo is again another solution to AG = and if u/q = u/q' < 1, 
Mq = M c , hence AG < 0. So, M depends on the ratio u/q, which is another indication 
that u and q are not totally independent. 

If J 7^ 0, G w still vanishes at M c , but G g vanishes at some other value M' different from 
M c . For j = and w/g > u/q', AG ^ for M ^ M and AG = for M = M . For j = 
and u/q < u/q' < 1, AG > for M > M > M min , where M min is the bound given by 
eq.(25) and q' in this case is given by AG = for M = M min = M ex , where M ex is the value 
of the extremal black hole such that M C 2 X = \(Q 2 + ^Q 2 + 4J 2 ). 

If J ^ and j ^ 0, if w/g > w/g', AG < for all M, while if u/q = u/q' < 1 and j < f, 
AG > for M > M > M min and AG < for M > M or M > M min . However, as j gets 
larger, the gap between M and M min shrinks, at that, if j > j' for some j', AG < for all 
M again. 

If J ^ and q = 0, for w/j > w/j', AG ^ for M ^ M , while for u/j < u/f, AG < 
for all M, where j' is determined by AG = for M = M ex . The G; = and J^0 case has 
a similar behavior, too. 

V. DISCUSSION 

We have shown the physical implication of the Davies critical phenomenon from the 
tunneling point of view. The sign of specific heat is related to the sign of the nonthermal 
contribution to the black hole radiation. We have also shown that the emission rate based 
on the tunneling peaks at (or near, if charge or angular momentum emissions are allowed) 
the Davies critical point. The argument based on the effective free energies show there are 
two phases, indicating whether the emission is spontaneous or not. 

Negative specific heat is common in collapsing self-gravitating systems which are isolated[14], 
so the real question to address in the black hole context is what happens to the black hole 
when the specific heat is positive. A black hole can never be truly a stable thermodynamic 
system because either matter always falls in and/or radiations always come out. The 



14 



equilibrium between falling matter and outward emission does not save the situation because 
of the case without matter available to fall and black hole solution does not distinguish it. If 
isolated, a black hole cannot evolve adiabatically without violating the charge conservation 
or angular momentum conservation. In case radiations come out at the same time matter 
falls, we cannot use tunneling argument to test the adiabatic evolution because the emission 
rate remains constant for any temperature, hence does not make sense. One could argue 
based on thermodynamic relations, using = dS = dM — V H dQ — Q H dJ, to check if 
there is a solution with a reasonable boundary condition. However, in the limit Q and J 
vanish, M has to be constant. Since M cannot be zero because, otherwise, it violates the 
unitarity, we end up with a Schwarzschild black hole whose specific heat is always negative. 
So the positive sign of the specific heat does not really imply its stability, unlike normal 
thermodynamic systems. Even in the extremal limit, it is not perfectly clear if the black hole 
is stable. (For example, see [15].) T H = 0, but we still have high charge-to-mass or angular- 
momentum-to-mass ratio so that it is not clear how it can sustain the zero temperature. 
We expect there should be no radiation, but the tunneling argument based on the WKB 
method fails in the extremal limit so that we cannot confirm, although it appears to be 
consistent. It will be interesting if one can check whether a black hole can really evolve into 
the extremal limit, or some kind of phase transition takes place at the extremal limit and it 
actually fails to be a black hole. 

From the tunneling point of view, the positive specific heat indicates that the radiation 
is suppressed by the nonthermal contribution. The effective Helmholtz free energy leads to 
AF > 0, indicating this suppression does not take place spontaneously. We suspect that this 
signals there should be another phenomenon taking place presumably right outside of the 
horizon since it may not be an intrinsic black hole phenomenon, which will end up increasing 
the black hole radiation toward the thermal level direction, like that of Damour-Ruffini[8] and 
Blandford-Znajek[9]. Therefore, we also suspect that there might be a connection between 
Parikh-Wilczek's nonthermal radiation[6] and that of Damour-Rufnni[8] and Blandford- 
Znajek[9]. The effort to relate them is in progress and will be reported elsewhere. 
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Appendix A: Yet Another Derivation of Hawking Temperature 

The black hole temperature can be obtained by performing the Wick rotation of the 
metric so that the period of the Euclidean time along a circle can be identified as the inverse 
temperature. Consider a generic metric with the Euclidean signature 

ds 2 = U(r)V(r)dt 2 + yy^dr 2 + rW. (Al) 

V(r) 

A godly traveler moving along a geodesic should be back to the original location after this 
time period, crossing the horizon back and forth. So we can compute the traveling period (3 

as 

'-jf*-/!*- 2 /,!*-'/,* ™ 

where P denotes a geodesic path crossing the horizon and we have eq.(ll) with t = t E - Note 
that for ds 2 = and 9 = = 

f = \F^- < A3 > 

dr V g tt 

So we might be tempted just to use this, but in the Euclidean signature this is imaginary 
so that we need to be a little bit more cautious and justify it. 

Without loss of generality and for convenience, we choose 9 = = 0, then nontrivial 
geodesic equations are 

= i+(\n(UV))'ir, (A4a) 
= f - \V{UV)'i 2 - \yf 2 , (A4b) 

where the dot denotes the derivative with respect to the affine parameter A and the prime 
denotes the derivative with respect to r. Eq.(A4a) can be integrated to 

UVi = ci, (A5) 

where c\ is an integration constant. Using this, eq.(A4b) can be written as 

°= f -^ v wr^-w f2 - < A6 > 
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Multiplying 2r/V, we obtain 



= 



d 

dX 



+ 4- 



d 



hence 



r 2 

'" 2 + V 



d\ \UV 



c 2 V, 



(AT) 



(A8) 



where c 2 is another integration constant. Using eq.(A5) and eq.(A8), we can perform the 
integration eq.(A2): 



, dt , 

= 2 d~r dr 

I p til 



dr 



Cl 



p UV y/czV-cj/U' 



(A9) 



Note that we can choose c 2 = for our geodesic path such that eq.(A3) is more or less 
justified. This integral is multi-valued, so we should choose the smallest nonvanishing value 
as the period. Then the integration simply becomes (assuming the pole is a simple pole, 
which is usually the case), using the Cauchy theorem after analytic continuation with a 
suitable path, 



(9 



i Jf 



dr—=— 

p Vuv 

2 r> 1 

-2ni Res 

l r=r H y/UV 



= Aii lim 



d_ 

dr 



Vuv 



(A10) 



The black hole temperature is, then, given by 



T, 



H 



1 



1 , d 
— hm — 

An r->r H dr 



— lim — 

An r->r H dr 




For U — 1, we have 



T, 



H 



1 



1 , d 

— hm —g tt . 

47T r-tr H dr 



(All) 



(A12) 



For KN black hole, set 9 = 0, which in turn sets g t( p = (also set = at the end), 
then the metric satisfies the same geodesic equations eqs.(A4a,A4b). So we can still apply 
eq.(All) to obtain the Hawking temperature for KN black hole as 



T, 



H 



r + — r_ 
An{r- 2 + + a 2 ) ' 



(A13) 
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Appendix B: Derivation of Eq.(20) for KN black hole 

In the WKB approximation, using the Hamilton- Jacobi theory, the emission of energy cu, 
charge q, and angular momentum j, is given by [11] [12] 

Im/ = Im J (h - p A A - dt, (Bl) 

where the generalized coordinates we use are (r, A, 0). Using eq.(18) and the fact that, as 
far as these generalized coordinates are concerned, S,Q, J are corresponding independent 
variables, we can obtain the following identities these generalized coordinates should satisfy: 



OH 1 OS ,„ , 

r = ^ = 7T^' B2a 

OPr PH OPr 

A = ™ = v* (B2b) 

Op A Op A 

■ dH „ dJ 

4> = — = n H — B2c 

optf, Opt, 

Substituting the integrations H — J dM, Pa = f dpA, and p<f, = f dp^, we obtain 

Im/ = Im J J y (rfM - VffdQ - fi^rfj) . (B3) 
In fact, using eq.(B2a), this is actually equivalent to 

Im/ = Im J p r rdt 

x f 9H , 
= Im / p r — — dt 

J OPr 

t ff, 1 dS , 
= Im / / dp r ——dt 
J J PHdp r 

^SS d ih iS - (B4) 

From eq.(Bl) and eq.(B4), we can consistently identify H as the Hamilton's characteristic 
function in this case. 

Finally, using eq.(ll), we can show that 



21m/ = AS = J d S = S{M - co,Q - q, J - j) - S(M, Q, J). 



(B5) 



We can assume J > without loss of generality and that jJ > 0. However, the formula 
is valid for any signs of Q and q, as long as AS < 0, hence we are not going to make any 
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assumption on the signs of Q and q at this moment. Let's expand AS in terms of u),q,j, 
such that 

dS dS .dS 
AS = - U dM- q dQ- 3 dJ 

&S_ 0^ 3*8 
+ 2 dM 2 + 2g dQ 2 2J dJ 2 

Q2g Q2g Q2g 

+ Uq dMdQ +q3 dQdl + 3U dJdM + (B6) 
Now define the following thermodynamic quantities for black holes: 

1 

ctQ = -p^.^f '■ charge expansion coefficient, (B7a) 
1 dJ 

a j = '■ angular momentum expansion coefficient, (B7b) 

kq = ^- '■ charge "compressibility", (B7c) 

k,j = — ( — J : angular momentum "compressibility". (B7d) 
J \oUh J t 

Then, from eq.(18), we can identify the following thermodynamic relations: 

' fa, (B8a) 
-PhV h , (B8b) 
-ftrftjf, (B8c) 



dM 
dS_ 

dQ 
dS 

dJ 

d 2 S 

dM 2 
d 2 S 
dQ 2 

dJ 2 

d 2 S d(3 H (3 2 H 



ft* 






Mc Q > 






^ M = 








ajJ 


Ph 



dMdQ dQ a Q Q 



dQdJ dQ" ' olqQ " S 

d 2 S _ d0 H fa 



(B8d) 
(B8e) 
(B8f) 
(B8g) 



d 2c ! d ft 2 ft 2 

^ (0 H n H ) = -^n H - ^-n H v H , (B8h) 



dJdM dJ ajJ 



(B8i) 
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Thus we have 



-2ImJ =AS 

= - (3 H (u - qV H - j^l h) 



1 2 (pH__fl L _ V \ 

2 Mc Q 2<1 \Qk q a Q Q H ) 
_ i . 2 ( Ph_ _ Ph n H 

2 \JKj OijJ 

- gj -tt^hVh — ll H 



-a/ — + h.o., (B9) 
J aj 



which is eq.(20). 



Appendix C: Derivation of Eq.(29) for KN black hole 

The effective Gibbs free energy is defined by 

G e fi = M eS — T eS S eS — V cS Q cS — Q eS J eff . (CI) 

M e ff and S e g can be identified as before in eq.(14), and the rest can be identified as following. 
The natural identification of AS in this case should be 

AS = -fa (w - qV cS - jft eff ) , (C2) 

then from eq.(B6) we can read off the effective quantities as 

_ dS 1 d 2 S ( d 2 S . d 2 S \ 

Pat ~dM~ 2U dM~ 2 ~ \ q dMdQ +3 dJdMj ' (C3a) 

n t, 9S 1 1 ■ 

PatVat = -dQ Jr ^W Jr ~ 23 dQdJ' (C3b) 

These identifications consistently truncate the effective Gibbs free energy to the effective 
Helmholtz free energy. And the definition of the effective Gibbs free energy eq.(Cl) also 
actually reduces to 

G cS = M- T cG S - V eS Q - Q eS J. (C4) 
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Now we can compute the difference between the effective Gibbs free energy and the usual 
Gibbs free energy at the same Hawking temperature as 

S 



where 



AG = G cS - G = -2- A/3 - QAV - JAQ, (C5) 



A/3 = & ff - f3 H 

= &(\-- + ~+ 3 --), (C6a) 
lH \ 2 Mc Q Qa Q Jaj> 

AV = V cS - V H 



! Q ( 1 ft, \ ! /Vff 1 

2 Q \kq a Q J \S a Q Q 



- ^A/3, (C6b) 
Ph 

AVL = Q cS - Q H 

l3 ( l far, \ , fl n /^H 1 



S a Q Q 



- ^A/3. (C6c) 
Ph 

This can be rewritten as 

AG = cjG u + qG q +jGj, (C7) 

where 



G ""-iiB^( 2M2+(32+ IS)' (C8a) 



where 



2M' 



(C8b) 



r - ^ fi i P3/2 2tt Q 2 + 2J 2 /M 2 / 
^ = (m5 - J~ H [ 2M +Q + B^))> (C8C) 



B = l — 2 , (C9a) 

A = 2B 3 / 2 + 2 - Q— A - 3^ + (C9b) 
M 4 M 2 M 6 V ; 
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